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摘 要 : 采集 塔克拉玛干 沙漠 北 缘 肖 塘 观 测 站 附近 的 土壤 ,应 用 Mastersizer -2000 激光 粒度 仪 测定 土壤 粒 径 ,按照 
伍 登 - 温 德 华 粒 级 标准 进行 土壤 成 分 分 类 ,发 现 土 壤 类 型 并 非 均 质 沙土 ,在 0 ~ 30 cm 深度 的 浅 层 存在 10% ~24% 


的 黏土 ,根据 土壤 成 分 比例 对 Noah 陆 面 模式 的 土壤 参数 重 计算 ,然后 利用 肖 塘 站 2011 年 3 月 22 日 至 7 月 25 日 的 


观测 数据 驱动 Noah 模式 ,对 比 修正 土壤 参数 前 后 的 模拟 效 呈 
量 偏差 减 小 10.47 W * m? , 感 热 通 量 偏差 减少 7.02 W - m”? ,地 表 净 辐射 偏差 减少 5.48 
W. m”, 且 效率 系数 达到 更 高 的 接受 程度 。 通 过 泰勒 图 分 析 了 土壤 参数 修正 对 Noah 模拟 效果 的 整体 作用 ,表明 


差 减少 0.62 C ,土壤 热 通 量 


Ro HERRI :更 新 土壤 参数 后 10 cm 土壤 温度 模拟 仿 


修正 土壤 参数 后 的 模式 中 热量 传输 扩散 有 一 定 的 优化 作用 ,10 cm 土壤 温度 和 热 通 量 土壤 温度 模拟 效果 得 到 改进 。 
关键 词 : 荒漠 ; 土壤 成 分 ; 陆 面 过 程 ; Noah 模式 ; 塔克拉玛干 沙漠 


陆 气 相互 作用 是 地 球 系统 中 对 陆 面 与 大 气 之 间 
的 物质 能 量 交 换 过 程 ,并 始终 影响 着 不 同时 空 尺度 
的 天 气 和 气候 变化 '"。 陆 面 为 天 气 和 气候 模式 提 
供 必 要 的 下 扑面 条件 ,大 气 强迫 陆 面 的 物理 特征 发 
生变 化 ,与 大 气 之 间 能 量 及 物质 的 交换 是 影响 天 气 
系统 的 重要 因素 ” 。 不 同类 型 下 热 面 的 陆 面 物理 
特征 有 明显 差异 ,因此 ,其 陆 面 过 程 也 存在 很 大 差 
异 "。 沙 漠 是 一 种 特殊 的 下 扑面 ,沙漠 地 表 反 照 率 
大 , 比 辐射 率 小 ,土壤 热 容量 小 ,含水 量 低 ,是 地 球 系 
统 的 重要 感 热源 ,对 区 域 能 量 平衡 及 气候 特征 具有 
重要 的 作用 ” 。 沙 漠 地 表 以 及 陆 面 过 程 与 全 球 其 
他 区 域 相 比 有 着 明显 的 差异 ,由 于 降水 量 极 低 , 植 被 
FAD , PETAL A PY LA ZAMAN TP, RPAH TE D UP HE 
区 地 气 间 能 量 交换 的 主要 方式 5 "。 在 土壤 热力 学 
研究 中 ,土壤 成 分 决定 了 其 始 .比热容 .热传导 率 、 热 
扩散 率 等 热力 学 性 质 “” 。 由 于 沙漠 地 表 的 独特 性 ， 
其 土壤 成 分 及 土壤 参数 是 影响 土壤 中 热量 传输 的 主 
要 因素 ,通用 的 陆 面 模式 中 沙漠 土壤 成 分 信息 和 参 
数 与 实际 值 存在 一 定 偏差 ,导致 陆 面 模式 对 沙漠 
下 垫 面 的 陆 面 过 程 模拟 效果 还 不 理想 "" 。 因 此 ,对 
沙漠 土壤 成 分 测定 和 陆 面 过 程 模 拟 分 析 是 一 种 提升 
沙漠 区 域 能 量 平 衡 模拟 效果 的 科学 方法 。 

近年 来 ,国内 外 开展 了 土壤 观测 和 陆 面 过 程 实 
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验 ,Liu 等 "基于 塔克拉玛干 沙漠 腹地 的 陆 气相 互 
作用 观测 试验 数据 ,修正 了 CoLM ( common land 
model) 模式 中 的 土壤 热传导 率 等 参数 ,能 量 通 量 模 
拟 效果 有 了 明显 提升 。 吴 小 波 等 "基于 Noah 陆 面 
过 程 模 型 模拟 青藏 高 原 植 被 和 土壤 特征 对 多 年 冻 土 
的 影响 , 发现 不 同时 期 的 土壤 状态 对 Noah 的 模拟 效 
果 有 明显 差异 。 土 壤 成 分 的 差异 会 导致 其 热力 学 参 
数 和 水 力学 参数 有 明显 差别 , 何 玉 洁 等 "对 青藏 高 
原 含 砂砾 石 土壤 导热 率 实 验 研 究 发 现 , 砂 砾石 对 土 
壤 导 热 率 有 显著 影响 ,提出 在 陆 面 过 程 模拟 研究 中 
必须 考虑 砂砾 石 对 土壤 热力 学 参数 的 作用 ,进而 提 
高 陆 面 过 程 模拟 的 精度 。 一 般 在 陆 面 模式 中 默认 土 
壤 为 均 质 分 布 ,会 导致 能 量 通 量 的 模拟 结果 有 不 确 
定性 。Han 4$? 3A Jy, CLM (community land mod- 
el) 的 土壤 参数 存在 不 确定 性 ,更 新 了 模式 中 的 土壤 
参数 后 ,CLM 模式 的 土壤 温 湿 度 、 热 通 量 等 的 模拟 
效果 有 一 定 改善 。 土 壤 性 质 在 不 同 的 空间 尺度 上 表 
现 出 高 度 的 非 均 质 性 ,其 准确 估算 混合 下 垫 面 的 土 
壤 参 数 仍然 是 一 个 重要 挑战 , Baroni 46") 提出 了 3 
种 描述 土壤 性 质 不 确定 性 的 微 扰 方 法 ,在 较 大 程度 
上 探索 出 土壤 参数 的 时 空 分 异 特征 ,进一步 完善 了 
水 文 模式 对 土壤 中 水 热 传 输 的 估算 。 曹 宾 琦 等 
对 塔克拉玛干 沙漠 北 缘 不 同 季 地 表 能 量 平 衡 闭合 特 
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征 进行 了 分 析 ,研究 表明 :地 表 以 下 5 cm 深度 的 能 
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火 青 等 :荒漠 土壤 成 分 对 陆 面 过 程 的 影响 


EEE 。 塔 克拉 玛 干 沙漠 也 是 中 国 极端 干旱 气候 


量 闭 合 程度 最 高 , 浅 层 土壤 的 能 量 通 量 对 陆 面 能 量 
平衡 非常 重要 。 不 同 土壤 成 分 的 热力 学 参数 对 陆 面 
模式 的 模拟 效果 有 重要 作用 ,这 些 参 数 既 能 影响 土 
壤 热 量 传 递 的 速度 和 深度 ,又 能 改变 土壤 水 分 传输 。 
沙漠 环境 不 便于 观测 实验 ,对 沙漠 区 域 土壤 成 分 及 
参数 的 测定 和 估算 存在 较 大 空白 ,本 项 目 利 用 Mas- 
tersizer — 2000 激光 粒度 仪 测定 塔克拉玛干 沙漠 北 
缘 肖 塘 观测 站 的 土壤 成 分 ,根据 成 分 比例 计算 Noah 
模式 层 土壤 参数 ,通过 Noah 模式 对 比试 验 分 析 , 探 
索 沙 漠 土壤 成 分 对 陆 面 过 程 的 影响 。 


1 研究 区 和 数据 


1.1 研究 区 概况 

塔克拉玛干 沙漠 是 中 国 面积 最 大 的 沙漠 (总 面 
积 3.38 x 10^ km? ,东西 相距 1 000 km、 南 北 相 距 
400 km) ,分 布 范围 37。~ 41°N、77° ~90°E, fir FS 
HUK x MEP RO” ,该 沙漠 大 多 为 流动 沙丘 和 沙丘 
链 , 如 300 m 高 的 金字 塔 形 沙丘 和 新 月 形 沙 


区 ,年 平均 降水 量 小 于 100 mm, 最 低 只 有 4 ~5 mm; 
而 年 平均 蒸发 量 达 3 400 mm, 有 显著 的 大 陆 性 气候 
REC? 。 

肖 塘 陆 气 相互 作用 观测 试验 站 (40°48'N， 
84*18'E, FIR 933 m, 简称 肖 塘 站 ) ,位 于 塔 克 拉 玛 
干 沙漠 北 缘 ( 图 1a) ,距离 塔里木 河 40 km, 位 于 平 
坦 的 风蚀 裸露 古河 床上 ,下 垫 面 为 自然 沙 面 ,周围 的 
新 月 形 沙 丘 及 沙丘 链 一 般 宽度 为 6~ 15 mm。 该 观测 
站 配置 有 100 m 的 11 层 梯 度 观 测 铁塔 \10 m 通 量 
观测 塔 7 层 土壤 温 湿度 和 通 量 观测 系统 ,以 及 9 种 
辐射 观测 仪 (图 1b)。 

1.0 数据 

于 2018 年 4 月 27 日 在 塔克拉玛干 沙漠 北 缘 新 
月 形 沙丘 及 平地 挖 断面 ,如 图 2a 所 示 。 新 月 形 沙丘 
及 平地 沙漠 南 缘 金字 塔 沙丘 及 绿洲 平地 采集 不 同 
部 位 不 同 深度 的 沙 粒 样本 开展 粒度 实验 分 析 。 其 
中 ,沙漠 北 缘 肖 塘 平 地 沙 粒 样本 采集 10 层 剖面 ,分 
Sig 0,3 .6 12 .20 ,40 .60 .100 .170 .200 cm; 其 余 各 


YE : 源 于 Google Earth, 
图 1 观测 站 位 置 示意 图 和 10 m 观测 塔 


Fig. 1 Locations of the observation station and the 10 — m-high observation tower 


2 土壤 成 分 采样 剖面 和 Mastersizer — 2000 激光 粒度 仪 


Fig.2 Sampling profile of soil composition and Masterizer — 2000 Laser Particle Sizer 
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沙丘 部 位 及 平地 均 和 采集 5 层 剖面 ,分 别 为 0.5、10、 
20,40 cm 深度 。 所 有 剖面 中 未 发 现 有 植被 根系 ,为 
HED ERFA LERENI 


2 研究 方法 


2.1 Noah 陆 面 模型 
Noah LSM 是 基于 OSU 陆 面 模型 发 展 而 来 , 遵 
循 强迫 -恢复 原理 。Noah 将 土壤 分 为 2 层 , 地 表 分 
为 单 层 积 雪 和 冠 层 。 基 于 热 扩散 方程 和 Richardson 
方程 分 别 计算 土壤 温度 和 土壤 水 含量 ,采用 Penman 
公式 计算 潜 热 通 量 ,其 植被 分 类 采用 USGS (U. S. 
geological survey ) 分 类 方案 ,土壤 数据 采用 世界 粮农 
组 织 全 球 土壤 数据 集 。 计 算 模 块 包括 :生物 气候 ,能 
量 、 冠 层 辐 射 传输 、 动 量 / 感 热 / 潜 热 通 量 交换 、 土 壤 
和 积 雪 温度 ,水文 、 动 态 植 被 等 模块 。Noah 能 够 全 
面 考虑 大 气 .植被 积 雪 等 局 部 因素 对 地 面 热 状态 的 
影响 ,描述 了 土壤 - 积 雪 - 植被 与 大 气 的 相互 过 程 ， 
能 够 模拟 土壤 温度 .土壤 含水 量 、 冠 层 含水 量 、 雪 深 、 
水 汽 、 能 量 通 量 . 向 上 长 短波 辐射 强度 等 。Noah 
模式 被 美国 环境 预报 中 心 广泛 应 用 于 天 气 /气候 模 
x. Haie zs WRF( weather research and fore- 
casting model) , RegCM ( regional climate model) 等 模 
式 ,经 过 应 用 和 检验 得 到 了 持续 的 发 展 ,主要 改进 
有 :土壤 分 层 改 为 4 层 (0.1.0.3.0. 6m 和 1.0 m); 
增加 了 和 冠 层 传导 方程 ””; 裸 地 蒸 散 和 植被 物候 
SEU ;新 的 径流 模式 ;热力 学 粗糙 度 的 参数 化 方 
案 及 地 表层 感 热 交换 系数 后 ) Noah (AE SE, EIR 
温度 计算 方程 可 以 表示 为 : 
Rig = (1 -a)S, *e(Ly -oT 4) (1) 


Raa =H+LE + G, (2) 


SUP RURSUM So 和 Ly 分 别 为 向 下 的 短波 、 长 
波 辐射 ,是 输入 数据 中 必要 的 驱动 要 素 ;7 是 地 表 
温度 ; o 是 Stefan-Boltzmann 常数 ;a 和 分 别 为 地 
表 反 照 率 和 地 表 比 辐射 率 。 在 能 量 平 衔 方 程 (2) 
中 ,五 是 感 热 通 量 ;LE 为 潜 热 通 量 ; Cu 为 土壤 热 通 
d. RPE H Æ Noah LSM 的 计算 可 以 表示 为 : 
H= - pe, Cu( Tair Ta) (3) 
式 中 ,为 空气 密度 ic, 是 空气 定 压 比 热 ;4 为 风速 ; 
7 为 距离 地 面 一 定 高 度 在 绝热 调整 后 的 空气 温度 ; 
7 为 地 表 温 度 ;C 为 感 热 交换 系数 ,是 一 个 非常 重 
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要 量 。 土 壤 热 通 量 Co 的 计算 方式 为 : 


T s. p T, 


Go =ky(O,) h (4) 


式 中 :fi 为 土壤 热传导 率 , 一 般 由 土壤 含水 量 (@) 
和 土壤 类 型 决定 ;7 为 土壤 上 层 温 度 ;h| 是 第 一 层 
土壤 的 中 间 深 度 。 在 干燥 的 土壤 条 件 下 ,@ 是 非常 
小 的 , 它 的 变化 在 大 多 时 候 基 本 可 以 忽略 ;kr TE il 
定 的 站 点 是 根据 土壤 参数 表 匹 配 的 固定 值 ,因此 ,在 
干旱 区 土壤 热 通 量 Cy 主要 依赖 于 hr 和 地 表 温 度 
7..。 土 壤 温 度 的 计算 如 下 : 


oT ð 
Ot | Oz 


式 中 :7 为 土壤 温度 ;i 为 时 间 ; 和 为 土壤 热传导 率 


[a Seco) | (5) 


(Wem + C7);C(9) 为 土壤 纹理 函数 ;z 为 土壤 
深度 。 


2.2 土壤 成 分 测定 方法 

沙 粒 样本 分 析 应 用 Mastersizer — 2000 激光 粒度 
仪 测定 (图 2b) ,测量 粒度 范围 为 0.02 ~ 2 000 pm。 
粒度 划分 标准 按照 伍 登 — 温 德 华 粒 级 标准 进行 , 粒 
度 参数 分 析 应 用 公式 (5 ) 计 算 平均 粒 径 。 

$= -logd (6) 

式 中 :d 是 沙 粒 直径 ;由 为 伍 登 - 温 德 华 -中 值 , 用 
于 对 土壤 质地 分 类 。 
2.3 ”分析 方 法 

在 本 次 研究 中 ,利用 2011 年 3 月 22 日 至 7 月 
24 日 的 观察 数据 驱动 Noah V3. 4 模式 ,这 期 间 肖 塘 
无 降水 ,有 多 云 / 沙 尘 等 天 气 ,为 客观 分 析 土 壤 成 分 
差异 对 陆 面 过 程 的 影响 , 需 晴好 天 气 , 因 此 ,选择 7 
月 1 一 10 日 模拟 结果 分 析 , 采 用 4 种 评 佑 指标 对 模 
型 的 模拟 结果 进行 分 析 , 分 别 为 :绝对 平均 误差 
(mean absolute error, MAE) ;效率 系数 (nash-sutcliffe 
efficiency, NSE) , NSE (EA fa 7625 #1) 1, 越 接近 1 表 
示 模 拟 结果 越 接近 观测 值 ,NSE = 1 时 ,表明 模拟 值 
与 观测 值 一 致 ,NSE > 0. 5 表示 模拟 效果 可 以 接受 ， 
当 NSE <0 时 则 表示 模拟 效果 差 ; 均 方 根 误差 ( root 
mean square error, RMSE) ,常用 于 模拟 值 与 观测 值 
之 差 的 平方 的 期 望 值 ;决定 性 系数 (及 ) ,用 于 表征 
模拟 值 与 观测 值 之 间 的 相关 性 。 其 中 ,前 3 个 指标 
主要 用 于 评价 模型 的 模拟 性 能 ,R* 用 于 泰勒 图 分 析 
模拟 效果 的 整体 性 能 。MAE、RMSE、NSE 和 R 分 
别 定义 为 : 


EKE SE reus ERRIA d Te D) Ug 


MAE =—¥ ( Ís, -0, ) (7) 
NSE =1 - $(5,-0,)?/ $(0,-0)" (8) 

RMSE = +$ (s, -0,)? (9) 
P= Benen] 

% (0; - 0)’ ESY (10) 


式 中 :5, 为 模拟 值 ;0, 为 观测 值 ; S 为 模拟 值 的 平均 
值 ; 0 为 观测 值 的 平均 值 ;n 为 样本 量 。 


3 ARNA 


3.1 土壤 成 分 测定 

在 剖面 采样 的 土壤 根据 伍 登 - 温 德 华 粒 级 标准 
粒度 划分 分 类 ,计算 出 每 一 层 的 土壤 成 分 比例 , 表 1 
为 多 层 土 壤 成 分 的 比例 ,经 统计 得 到 模式 层 的 土壤 
成 分 比例 (图 3) ,其 中 第 一 层 砂 质 占 75.47% , Bi 
24.53% ,第 二 层 砂 质 占 87. 8% ,黏土 12. 2% ,下 两 
层 均 为 砂 质 土壤 ,然后 按照 Noah 模式 中 沙 地 和 黏土 
混合 层 的 参数 , 见 表 2。 通 过 表 2 可 以 了 解 到 ,重新 
测定 的 土壤 参数 有 和 较 大 差异 ,土壤 热力 学 参数 减 小 ， 
土壤 水 力学 参数 增 大 。 
3.2 土壤 温度 

土壤 温度 (5S7) 是 陆 面 过 程 重 要 的 物理 量 ,直接 
影响 着 土壤 中 水 分 \ 水 汽 的 保持 和 运 移 ,对 土壤 中 的 
物理 过 程 均 有 有 影响。 图 4 给 出 了 Noah 的 10 cm + 
二 温度 观测 值 与 2 组 模拟 对 比 曲线 , Noah 模式 均 能 
够 较 好 地 反映 塔 中 站 的 土壤 温度 变化 ,原始 土壤 参 
数 模拟 的 土壤 温度 早 正 午时 明显 偏 高 ,在 夜间 有 偏 
低 现 象 ,修正 土壤 参数 后 的 模拟 结果 较 好 地 改善 了 
这 种 情况 。 根 据 表 3 可 知 ,修正 后 的 土壤 热传导 率 
系数 明显 减 小 ,其 绝对 平均 偏差 从 1. 332 CWE 
0.712 °C ,RMSE 由 1.56 %C 减 小 到 0.78 °C ,模拟 效 


率 系数 也 有 明显 提升 。 土 壤 的 热量 来 源 是 太阳 辆 
射 ,每 种 土壤 吸收 一 定 热量 后 温度 变化 的 程度 存在 
差异 ,主要 是 其 土壤 热力 学 参数 不 同 ,通过 表 2 可 看 
出 , 按 砂 质 和 黏土 比例 计算 的 热传导 系数 减 小 了 
0.36 ,根据 土壤 温度 计算 公式 (5) 可 知 ,土壤 热传导 
率 的 减 小 会 导致 土壤 温度 振幅 减 小 ,改善 了 中 午 高 
估 和 夜间 的 低估 程度 。 
3.3 土壤 热 通 量 

土壤 热 通 量 Gy 作为 地 表 能 量 平衡 方程 中 的 重 
要 分 量 之 一 ,是 表征 地 表 接 收 到 太阳 辐射 癌 下 传 到 
下 层 土壤 的 热量 ,一般 定义 为 向 下 传导 热量 为 正 通 
量 ,反正 为 负 。 土 壤 热 通 量 的 计算 主要 取决 于 地 表 
温度 ,土壤 热力 学 参数 .土壤 水 含量 。 通 过 图 5 可 以 
看 出 , Noah 能 够 较 好 模拟 出 土壤 热 通 量 的 变化 趋 
势 , 但 原始 参数 模拟 值 在 中 午 明显 偏 高 ,最 高 超过 观 
测 200W + m”。 根 据 土壤 成 分 重新 计算 土壤 参数 


表 1 多 层 土壤 成 分 比例 


Tab.1 composition ratios of multi-layer soil 
粗 砂 Pe aque» — 粗 粉 砂 ” 细 粉 砂 Rit 
0 cm 0 0.55 17.79 12.45 47.72 21.49 
3 cm 0 0.11 0.98 12.52 50.78 35.62 
6 cm 0 0 1.74 28.72 41.98 27.55 
12 cm 0 0 2.45 23.7 44.33 29.51 
20 cm 0 0.69 56 14.45 16.66 12.2 
40 cm 0 4.22 93.9 1.87 0 0 
60 cm 0 4.02 94.09 1.89 0 0 
100cm 0 5.69 92.4 1.91 0 0 
平均 值 0 2.425 53.88 10.909 20.147 12.637 


0.1-0.8 m Rb £87.80% 黏土 12.20% 


0.3~0.6 m 砂 土 100% 


0.6~1.0 m 砂 土 100% 


图 3 土壤 分 层 及 成 分 信息 示意 图 


Information of soil stratification and composition 


Fig.3 


表 2 土壤 参数 
Tab.2 Soil parameters 


Sand BB DRYSMC F11 MAXSMC REFSMC SATPSI WLTSMC QTZ 
SP-Org 2.79 0.011 -0.472 0.339 0.236 0. 069 0.01 0.92 
SP-Res 4.67 0.037 — 0.832 0.339 0.274 0.154 0.037 0.77 


TE: SP-Org 为 原始 土壤 参数 ;SP-Res 为 重新 测定 后 土壤 参数 。BB 为 计算 水 势 的 B 参数 ;DRYSMC Oh E E BEAR Ae BUB; FLL 为 土壤 热传导 系数 ; 


MAXSMC 为 饱和 土壤 湿度 ;REFSMC 为 水 分 传输 时 相对 土壤 湿度 ;SATPSI 为 土壤 基质 势 伯 和 湿度 ; WLTSMC 为 冰冻 层 熔 融 点 湿度 ;QTZ 为 土壤 


石英 含量 。 
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To" 
—— ST-Obs( 观 测 ) 一 一 ST-Org( 原 始 ) 
p 45 -一 ST-Res( 修 正 土壤 参数 后 模拟 ) 
E 
mH 
ER 
Hc 
g 
2 
S 
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日 期 /月 -日 
到 4 10 cm 土壤 温度 模拟 分 析 曲 线 
Fig.4 Simulated result of soil temperature at 10 — cm depth 


表 3 误差 统计 分 析 


Tab.3 Statistical results of errors 


指标 / MAE RMSE NSE 
变量 Org Res Org Res Org Res 
ST 1.332 0.712 1.56 0.78 0.72 0.93 
Go 39.28 28.81 61.63 45.88 0.47 0.71 
H 18.67 11.65 9.32 9.30 0.98 — 0.99 
Raa 28.83 23.35 31.40 15.09 0.97 0.99 


TE: Org 为 原始 土壤 参数 模拟 ; Res 为 修正 土壤 参数 后 模拟 。 


Gi-Obs( 观 测 ) Gu-Org( 原 始 ) 


个 Gi-Res( 修 正 土壤 参数 后 模拟 ) 
= 400 
= 300 
S 200 
til 109 } 
时 0 
ER -100 
4H = Nn o + © o F o Q o 
o à &©& à ð à S à & J 
iiiiiiiiis 
日 期 /月 -日 
ALS 土壤 热 通 量 模拟 分 析 曲 线 图 
Fig.5 Simulated result of soil heat flux 


后 的 模拟 结果 有 一 定 改善 ,绝对 平均 偏差 减 小 
10.47 W * m? ,RMSE 减 小 15.75 W * m ? ,但 仍然 
高 估 。 考 虑 到 肖 塘 站 长 期 没有 降水 ,土壤 极其 干燥 ， 
由 于 石英 的 比热容 小 且 热 传导 率 较 高 ,修正 后 土壤 
参数 石英 (QZT) 成 分 减少 ,所 以 ,表层 混合 土壤 计算 
的 热传导 系数 减 小 。 根 据 土壤 热 通 量 计算 公 式 (4) 
可 知 ,导致 土壤 热 通 量 高 估 减 小 的 原因 主要 是 土壤 
热力 学 参数 的 改变 , 减 小 了 热量 向 下 层 土壤 传递 速 
度 与 热量 。 
3.4 miim 

感 热 通 量 有 是 地 表 能 量 传输 的 主要 分 量 , 是 塔 
克拉 玛 干 沙漠 的 主要 热源 ,能 反映 出 地 面 与 大 气 之 
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间 的 热量 交换 强度 ,是 评价 陆 面 过 程 模式 模拟 效果 
的 关键 物理 量 。 从 图 6 可 知 , 感 热 通 量 模 拟 值 均 能 
够 较 好 地 反映 其 变化 趋势 ,原始 土壤 参数 的 感 热 通 
量 模 拟 值 在 中 午 有 局 部 偏 高 ,修正 土壤 参数 后 有 微 
弱 改 善 ,绝对 平均 偏差 减 小 了 7. 02 Wem? YAY 
根 误差 和 效率 系数 提升 效果 甚 微 , 均 没 有 反映 出 7 
AS 日 夜间 感 热 通 量 的 剧烈 波动 。 测 定 土壤 参数 可 
知 ,土壤 热传导 系数 有 所 减 小 ,根据 感 热 通 量 计算 公 
式 (3) ,沙漠 区 域 的 感 热 通 量 主要 取决 于 地 表 温 度 
和 土壤 温度 差 , 由 上 已 知 , 土 壤 温度 的 高 估 得 到 改 
善 , 因 此 ,在 正午 时 的 感 热 通 量 高 估 有 所 降低 。 


一 一 H-Obs( 观 测 ) 


一 一 HH-Org( 原 始 ) 


E 400 一 一 H-Res( 修 正 土壤 参数 后 模拟 ) 

(= 

+ 300 

E 

So 

到 100 

= 0 
一 | e [o^ ~ 19 © bt oo 2 © 
TT 
也 i Lr i p r- © © © © 
o 一 e e =) © e © © e 


日 期 /月 -日 
图 6 感 热 通 量 模 拟 分 析 曲 线 图 


Fig.6 Simulated result of sensible heat flux 


3.5 #48 

地 表 净 辐射 R,, 是 地 表 吸 收 的 太阳 辐射 能 量 ， 
是 估算 地 表 能 量 平衡 的 关键 要 素 ,对 天 气 / 气 候 预 测 
有 重要 的 作用 ,也 是 评估 太阳 能 资源 的 重要 指标 。 
图 7 为 Noah 净 辐 射 模拟 值 与 观测 值 对比 曲 线 图 ,2 
组 净 辐 射 模拟 值 的 变化 趋势 与 观测 值 非常 接近 ,但 
是 原始 土壤 参数 的 模拟 结果 在 正午 有 偏 高 ,平均 偏 
高 50 W- m ,修正 土壤 参数 后 的 模拟 结果 改善 了 
这 种 情况 ,其 模拟 结果 基本 和 观测 值 吻 合 。 从 误差 
统计 可 知 ,修正 土壤 参数 后 净 辐 射 模拟 结果 的 绝对 
平均 偏差 减 小 了 5.48W . m ^ , RMSEJR/]N f 16.31 


一 R,,,-Obs(Miill) — R,,-Org Rt) 


EM — R -Res( 修 正 土 壤 参 数 后 模拟 ) 

=} 

Bs 

dR 

m 

ab 
一 N [n^] bu [te] Ke] œ~ oo © © 
© © © © © © c c So = 
ee a ae ae 
© © © © © © © © © © 
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图 7 净 辐 射 模拟 分 析 曲 线 图 


Fig.7 Simulated result of net radiation 
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W- m”, OR AB NSE 达到 更 高 的 接受 程度 。 干 
燥 环 境 的 潜 热 通 量 基本 可 以 忽略 ,通过 能 量 平衡 公 
式 (2) 可 知 , 净 辐 射 结 果 取 决 于 感 热 通 量 和 土壤 热 
通 量 的 计算 ,因此 ,符合 环境 的 土壤 参数 ,有 利于 能 
量 的 准确 估算 。 
3.6 小 结 

综合 以 上 对 土壤 温度 .土壤 热 通 量 、 感 热 通 量 、 
净 辐 射 的 2 组 模拟 结果 分 析 , 本 次 研究 加 入 泰勒 图 
对 模型 模拟 性 能 进行 了 评估 ,泰勒 图 是 一 种 能 从 多 
方面 评估 模型 模拟 性 能 的 分 析 图 ,可 以 同时 反映 相 
关 性 和 误差 情况 。 图 8 为 土壤 参数 修正 前 后 Noah 
的 模拟 性 能 评估 泰勒 图 , 归 一 化 标准 差 用 直径 表示 ， 
相关 性 用 角度 表示 , 当 点 越 靠 近 横 轴 , 则 相关 性 越 
高 。 从 图 中 可 以 看 出 ,土壤 温度 改进 效果 最 明显 , 感 
热 通 量 提升 效果 甚 微 ,土壤 热 通 量 模拟 结果 的 相关 
性 有 降低 。 但 是 总 体 来 看 ,模式 对 能 量 通 量 的 模拟 
效果 有 改善 ,说 明 符合 环境 的 土壤 成 分 信息 有 助 于 
提高 模式 模拟 性 能 。 


St ow 
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图 8 泰勒 图 


Fig.8 Taylor diagram 


4 结论 


在 陆 面 模式 中 土壤 参数 一 般 由 参数 表 匹 配 , 认 
为 是 均值 类 型 的 土壤 ,然而 在 实际 中 存在 混合 土壤 
类 型 情况 较 多 ,本文 研究 对 肖 塘 的 土壤 成 分 采集 和 
粒 径 分 类 ,测定 不 同 层 的 土壤 成 分 和 比例 ,根据 沙 和 
黏土 成 分 比例 重新 计算 了 Noah 的 土壤 参数 ,并 利用 
该 站 的 观测 数据 驱动 作 了 2 组 陆 面 过 程 模拟 和 分 
Br ,得 到 以 下 几 个 结论 : 


(1) 塔克拉玛干 沙漠 北 缘 沙漠 土壤 极其 干燥 ， 
但 其 中 不 仅 只 有 沙 质 ,在 0 ~30 cm 深度 的 浅 层 存 在 
10% ~24% AREE , A ERIS FE BA , ARE ABH E 
成 分 比例 重新 计算 土壤 参数 后 发 现 ,混合 土壤 的 参 
数 有 较 大 差异 ,其 中 表层 土壤 热传导 率 减 小 ,其 他 土 
壤 水 力学 参数 增 大 。 

(2) 在 更 新 土壤 参数 表 后 ,Noah 模式 对 10 cm 
土壤 温度 的 模拟 效果 有 明显 提升 ,修正 后 的 土壤 热 
传导 率 系数 减 小 ,改善 了 土壤 温度 在 中 午 和 夜间 偏 
差 较 大 现象 ,在 这 期 间 无 降水 ,土壤 水 分 基本 没 变 
化 ,土壤 水 力学 参数 的 增 大 不 会 对 热 传 输 过 程 产生 
影响 , 谨 明 模式 默认 的 沙 质 土壤 的 热传导 率 偏 高 , 修 
正 后 的 土壤 热传导 率 减 小 ,能 够 减缓 热量 向 下 传递 
的 速度 ,因此 , 减 小 了 温度 变化 的 振幅 。 

(3) 更 新 土壤 参数 表 后 ,表层 土壤 热传导 系数 
有 所 减 小 ,土壤 热 通 量 在 修正 土壤 参数 后 的 模拟 结 
果 一 定 程 度 上 减 小 了 中 午 模拟 偏 高 的 现象 ,表明 土 
培 成 分 对 热量 向 下 传输 有 明显 作用 。 正 午时 的 感 热 
通 量 高 估 有 所 降低 ,有 微弱 改善 ,但 均 方 根 误差 和 效 
率 系数 提升 效果 其 微 ,说 明 地 表 感 热 通 量 对 土壤 热 
力学 参数 不 是 非常 敏感 。 另 外 ,修正 土壤 参数 后 , 净 
辐射 模拟 结果 与 观测 值 更 加 吻合 , 净 辐 射 的 模拟 效 
果 提 升 取决 于 感 热 通 量 和 土壤 热 通 量 的 结果 。 误 差 
统计 指标 表明 均 有 改进 ,因此 ,修正 土壤 参数 更 符合 
肖 塘 的 环境 ,有 利于 提高 地 表 能 量 平衡 闭合 率 。 

综 上 可 知 , 模 式 上 默认 每 一 类 土壤 为 均 质 成 分 , 真 
实 环境 中 表层 土壤 一 般 会 混合 多 种 土壤 成 分 , Noah 
模式 的 土壤 参数 表 并 不 符合 实际 情况 ,本 研究 表明 ， 
更 符合 环境 特征 的 土壤 成 分 信息 有 助 于 提高 模式 对 
热 传 输 模 拟 性 能 。 但 是 大 面积 区 域 土壤 属性 的 时 空 
异 质 性 测定 通过 剖面 采样 观测 是 不 现实 的 ,探索 多 
光谱 遥感 和 地 面 观 测 联合 反 演 多 种 土壤 参数 是 未 来 
面向 区 域 研究 和 应 用 的 科学 方法 。 


参考 文献 ( References ) : 

(1) Dickinson R E. Land-atmosphere interaction J). Reviews of Geo- 
physics, 1995 ,33 (33) :917 - 922. 

(2) Chen F,Zhang Y. On the coupling strength between the land surface 
and the atmosphere ; From viewpoint of surface exchange coefficients 
(JJ. Geophysical Research Letters ,2009 ,36 (36) :207 -213. 

(3) Wei J, Dirmeyer P, Guo Z ,et al. How much do different land mod- 
els matter for climate simulation? Part I; Climatology and Variabili- 
ty CJ). Journal of Climate ,2010 ,23 (11) :3 135 -3 145. 

(4) Yang T, Wang X, Zhao C, et al. Changes of climate extremes in a 


202001.00040v1 


chinaXiv 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


TOS 区 


typical arid zone; Observations and multimodel ensemble projec- 
tions( J). Journal of Geophysical Research Atmospheres ,2011 ,116 
(D19) :106 - 124. 
Rodwell M J, Hoskins B J. Monsoons and the dynamics of deserts 
(JJ. Quarterly Journal of the Royal Meteorological Society , 1996 , 
122(534) :1 385 - 1 404. 

1, 李 振 朝 ,等 . 西北 干旱 区 戈壁 下 垫 面 陆 面 过 程 特 
HERJER IUDICI]. 太阳 能 学 报 ,2016 ,37(9) :2 422 - 2 429. 
( Wang Chao, Wei Zhigang, Li Zhenchao,et al. Simulation analysis 


pH + 
-EJEB , TRR 


of land surface process over Gobi at arid region Northwest China 
during summer( J]. Acta Energiae Solaris Sinica, 2016, 37 (9): 
2 422 -2 429.) 

ZEA, THES, (EEO, SE. 金 塔 绿洲 的 辐射 平衡 特征 和 地 表 
FEET ESE (I). 太阳 能 学 报 ,2009,30(12):1 614 - 1 620. (Li 


Wanli, Lii Shihua, Fu Shengming, et al. A study on radiation bal- 


ance and surface energy in Jinta Oasis( J]. Acta Energiae Solaris 
Sinica ,2009 ,30(12) :1 614 - 1 620. J 

Alkhaier F, Flerchinger G N,Su Z. Shallow groundwater effect on 
land surface temperature and surface energy balance under bare 
soil conditions ; Modeling and description( J}. Hydrology and Earth 
System Sciences ,2012 ,16(7) :1 817 -1 831. 

李 火 青 , 吴 新 萍 , 买 买 提 艾 力 . 买 买 提 依 明 ,等 . 塔克拉玛干 沙 
漠 地 表 浅 层 土壤 热 扩 散 率 温度 和 热 通 量 计算 方法 的 比较 研 
究 [J]. 土壤 通报 ,2016,47(4):805 -813. ( Li Huoqing, Wu 


Xinping, Ali Mamtimin, et al. Comparation between the different 
methods toestimate soil thermal diffusivity, temperature and heat 
Flux in Sub-surface of Taklimakan Desert( J]. Chinese Journal of 
Soil Science ,2016 ,47(4) :805 - 813. ) 

郑 辉 , 刘 树 华 . vb it Bali TM ob Ee UG RA CJ). 地 球 物理 学 
报 ,2013 ,56(7) :2 207 -2 217. ( Zheng Hui, Liu Shuhua. Land 


surface parameterization and modeling over desert [ J]. Chinese 
Journal of Geophysics ,2013 ,56(7) :2 207 -2 217. J 

Liu Y Q, He Q, Ali Mamtimin. Improving the CoLM in Taklimakan 
desert hinterland with accurate key parameters and an appropriate 


parameterization scheme( J}. Advances in Atmospheric Sciences, 


2012 ,29 (2) :381 — 390. 
吴 小 波 , 南 卓 铜 , 王 维 真 ,等 . 基于 Noah 陆 面 过 程 模型 模拟 青 
藏 高 原 植被 和 土壤 特征 对 多 年 冻 土 的 影响 [ 相 .: 冰川 冻 土 ， 
2018,40 (2 ) : 280 - 286. ( Wu Xiaobo, Nan Zhuotong, Wang 


Weizhen ,et al. Simulation of the impact of vegetation and soil char- 
acteristics on permafrost over the Tibetan Plateau based on the No- 
ah land surface model( J). Journal of Glaciology and Geocryology, 
2018 ,40(2) :280 - 286. J 

fu] dci EAE, SHEER. 青藏 高 原 含 砂砾 石 土壤 导热 率 实 验 
WE CI). 冰川 冻 土 ,2017,39(2):119 - 126. ( He Yujie, Yi 


Shuhua, Guo Xinlei. Experimental study on thermal conductivity of 
soil with gravel on the Qinghai—Tibet Plateau( J}. Journal of Glac- 
iology and Geocryology ,2017 ,39(2) :119 - 126. ) 


(14) 


(15) 


(16) 


[17] 


(18) 


(19) 


(20) 


(21) 


(22) 


(23) 


(24) 


[25] 


ChinaXiv 合 作 期 刊 


Han X, Franssen H, Montzka C ,et al. Soil moisture and soil proper- 
ties estimation in the Community Land Model with synthetic bright- 
ness temperature observations( J]. Water Resources Research ,2014 , 
50(7) :6 081 -6 105. 

Baroni G, Zink M, Kumar R, et al. Effects of uncertainty in soil 
properties on simulated hydrological states and fluxes at different 
spatio-temporal scales {J}. Hydrology & Earth System Sciences, 
2017 ,21(5) :1 -37. 

DE Hay AMT, aA T. 塔克拉玛干 沙漠 北 缘 夏秋 冬季 地 表 
能 量 平衡 闭合 特征 [机 . 干旱 区 研究 ,2018,35 (4) :830 - 839. 
( Cao Huanqi, He Qing, Jin Lili, et al. Surface energy balance clo- 


sure of summer, autumn and winter in the Northern Margin of Takli- 
makan Desert( J). Arid Zone Research ,2018 ,35 (4) :830 — 839. ] 
Sun J. The age of the Taklimakan Desert( J). Science ,2006 ,312 
(5 780) :1 621 - 1 621. 

钱 正安 , 宋 敏 红 , 李 万 元 . 近 50 年 来 中 国 北方 沙尘暴 的 分 布 及 
BE (boy WT CI). 中 国 沙漠 ,2002,22(2):106 - 111. ( Qian 
Zhengan , Song Minhong, Li Wanyuan, et al. Analysis of distribution 
and trend of sandstorm in North China in recent 50 years( J}. Jour- 
nal of Desert Research ,2002 ,22(2) :106 - 111. J 

魏 文 寿 . 沙漠 气象 学 [ M]. 北京 :气象 出 版 社 ,2008 : 197. ( Wei 
Wenshou. Desert Meteorology ( M ). Beijing: Meteorological Publish- 
ing House ,2008 : 197. ] 

金 莉莉 , 李 振 杰 , 何 清 ,等 . CoLM 模式 对 塔克拉玛干 沙漠 北 缘 
陆 面 过 程 模拟 评估 及 修正 [站]. 中 国 沙漠 ,2018 ,38(6) :173 - 
188. (Jin Lili, Li Zhenjie, He Qing, et al. Evaluation on modeling 


of land-atmosphere interaction in the Northern marginal zone of 
Taklimakan Desert with CoLM and correction[ J}. Journal of Desert 
Research ,2018 ,38 (6) :173 —188. ] 

Chen F, Janji € Z, Mitchell K. Impact of atmospheric surface-layer 
parameterizations in the new land-surface scheme of the NCEP me- 
soscale Eta model( J). Boundary-Layer Meteorology ,1997 ,85(3) : 
39] - 421. 

Chen F, Mitchell K, Schaake J, et al. Modeling of land surface e- 
vaporation by four schemes and comparison with FIFE observations 
(JJ. Journal of Geophysical Research Atmospheres, 1996, 101 
(D3) :7 251 -7 268. 

Betts A K. Assessment of the land surface and boundary layer mod- 
els in two operational versions of the NCEP Eta model using FIFE 
data( J]. Monthly Weather Review, 1997, 125 (11) :2 896 - 
2 916. 

Schaake J C, Koren V , Duan Q Y ,et al. Simple water balance mod- 
el for estimating runoff at different spatial and temporal scales( J). 
Journal of Geophysical Research Atmospheres, 1996 , 101 ( D3 ) : 
7 461 -7 475. 

Chen Y Y, Yang K,He J, et al. Improving land surface temperature 
modeling for dry land of China( J). Journal of Geophysical Research 
Atmospheres ,2011. https: //doi. org/10. 1029/2011JD015921. 


EKE SE M ERRIA Bd V E Fe D) Ug 


Effects of Desert Soil Composition on Land Surface Processes 


LI Huo-qing , JIN Li-li', ALI Mamtimin', ZHAO Jia-wei 
(1. Institute of Desert Meteorology , China Meteorological Administration , Urumqi 830002 , Xinjiang , China ; 
2. Xinjiang Normal University , Urumqi 830013 , Xinjiang , China ) 


Abstract; In this study, the values of soil composition at Xiaotang Observation Station in the northern margin of 
the Taklimakan Desert were measured by Mastersizer — 2000 Laser Particle Sizer and analyzed. Soil composition was 
classified according to the Wooden-Wendward Particle Size Standard. It was found that the soil was not the homoge- 
neous sandy soil, and there was 10% -24% clay in the shallow soil layer at 0 - 30 — cm depth. The soil parameters 
in the Noah Land Surface Model were recalculated according to the proportions of soil compositions. The data ob- 
served at Xiaotang Station from March 22 to July 25,2011 were used to drive the Noah model ,and the simulated re- 
sults before and after amending the soil parameters were compared. The results showed that , after updating the soil 
parameters, the simulation effect of soil temperature was improved, the deviations of soil temperature , soil heat flux, 
sensible heat flux and surface net radiation were reduced by 0.62 °C ,10.47 W + m ?,7.02 W - m ^ and 5.48 W 

m ^ respectively , and the efficiency coefficient was more acceptable. The overall effect of soil parameter modifica- 
tion for the Noah simulation was analyzed by Taylor diagram. The results showed that the heat transfer and diffusion 
in the modified soil parameter model could be improved to a certain extent, and the simulation effect of soil tempera- 
ture could be improved obviously. 


Key words: desert; soil composition; land process; Noah model; Taklimakan Desert 


